Non-trivial Vacuum Effects in the 2-d 0(3) Model* 



M. Bogli", F. Niedermayer", M. Pepe'', and U.-J. Wiese" 
" Albert Emstein Center for Fundamental Physics, Institute for Theoretical Physics, 
Bern University, Sidlerstr. 5, 3012 Bern, Switzerland 
INFN, Istituto Nazionale di Fisica Nucleare, Sezione di Milano-Bicocca 

Edificio U2, Piazza della Scienza 3 - 20126 Milano, Italy 
* Dedicated to Peter Hasenfratz on the occasion of his 65th birthday 

We study S-vacua in the 2-d lattice 0(3) model using the standard action and an optimized 
constraint action with very small cut-off effects, combined with the geometric topological charge. 
Remarkably, dislocation lattice artifacts do not spoil the non-trivial continuum limit a.t 6 ^ Q, and 
there are different continuum theories for each value < 6* < tt. A very precise Monte Carlo study 
of the step scaling function indirectly confirms the exact S-matrix of the 2-d 0(3) model at ^ = tt. 



As Bethe showed in 1931, the spin i antiferroniag- 
netic Heisenberg chain is gapless [l|. In 1983 Haldane 
conjectured that the spin S chain has a gap for inte- 
ger S — 1,2,3,..., but is gapless for half-integer spins 
5' = ^,|,|,...[2|. In the semi-classical large S limit, he 
showed that the corresponding low-energy effective field 
theory is the 2-d 0(3) model at vacuum angle = 2ttS, 
with coupling = 2/ S. Chains with integer spin hence 
correspond to 6 — Q, while half- integer spin chains have 
9 = TT. The 2-d 0(3) model dX 9 = n should reduce to 
the k = \ WZNW model at low energies. Haldane's 
conjecture has been confirmed by numerical simulations 
of integer and half-integer quantum spin chains [6, 7]. A 
direct comparison with the 2-d 0(3) model at 6* = tt is 
more difficult, due to the notorious complex action prob- 
lem in numerical simulations. However, using a meron- 
cluster algorithm (an extension of the Wolff cluster algo- 
rithm [^) it was possible to simulate at 9 = n, and consis- 
tency with the WZNW model predictions was obtained 
within statistical errors In this paper, we also use a 
variant of a method developed by Hasenbusch l^, U | to 
simulate 0- vacuum effects in the 2-d 0(3) model with un- 
precedented per mille level precision. For the first time, 
this numerically confirms the conjectured exact S-matrix 
of the 2-d 0(3) model at = tt [12| beyond any rea- 
sonable doubt, which also implies that the model indeed 
reduces to the WZNW model at low energies. 

The 2-d 0(3) model is e quiv alent to the CP(1) model. 
2-d CP(iV - 1) models [H,] |14| share many features with 
4-d non-Abelian Yang-Mills theories: they are asymptoti- 
cally free, have a non-perturbatively generated massgap, 
instantons, as well as non-trivial 6'-vacua. CP{N — 1) 
models on the lattice have been studied by Berg and 
Liischer 15| who introduced a geometric definition for 
the lattice topological charge Q. Field configurations 
with Q = 1 and a minimal value of the lattice action 
are known as dislocations. In general, dislocations have 
a smaller action than Q = 1 instantons — the minimal 
action configurations of the continuum theory. When the 
dislocation action is less than a critical value (determined 
by the 1-loop /3-function coefficient), a semi-classical (but 
non-rigorous) argument suggests that the topological sus- 



ceptibility xt — {Q'^)/V (where V is the space-time vol- 
ume) should suffer from an ultra-violet power-law diver- 
gence in the continuum limit [16]. In CP{N — 1) models 
with A'^ > 4, the dislocation problem does not arise when 
one uses the standard lattice action in combination with 
the geometric definition of the lattice topological charge. 
In the CP(2) model (i.e. for = 3) the problem is ex- 
pected to arise, but can be avoided by the use of an im- 
proved lattice action [l7[, which pushes the dislocation 
action above the critical value. Finally, in the CP(1) 
(i.e. the 0(3)) model the critical value of the disloca- 
tion action agrees with the continuum instanton action. 
Consequently, dislocations can be suppressed only by a 
delicate fine-tuning of the lattice action. This was real- 
ized in a study with a classically perfect lattice action, 
which indeed led to a logarithmic rather than a power- 
law divergence of Xt El8|. Semiclassically, the logarithmic 
ultra-violet divergence also arises directly in the contin- 
uum [igI . [igj and may thus not be a lattice artifact. In 
that case, Xt would not be a physically meaningful quan- 
tity in the 2-d 0(3) model. Along the way, the whole 
concept of distinct topological sectors and corresponding 
6'-vacua has been questioned in the 2-d 0(3) model. In 
fact, one may suspect that 9 is an irrelevant parameter 
that renormalizes to zero non-perturbatively. In this pa- 
per, for the first time we demonstrate with high accuracy, 
that 9 is actually relevant and that each value < 9 < ir 
is associated with a different continuum theory. 

A first indication that dislocations may not have a dev- 
astating effect on the continuum limit arose in a recent 
study of the 2-d 0(3) model at 6* = with a topologi- 
cal lattice action [20|, where we found again just a loga- 
rithmic divergence of Xt ■ Topological lattice actions are 
invariant against small deformations of the lattice fields. 
In particular, in [20| an action with a constraint on the 
maximally allowed angle between neighboring spins has 
been used [2l[. All allowed configurations (i.e. those that 
respect the constraint) have the same action value zero. 
As a consequence, this lattice model does not have the 
correct classical continuum limit, it violates the Schwarz 
inequality between action and topological charge, and it 
cannot be treated in perturbation theory. Despite these 
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FIG. 1. Triangulated square lattice: the triangles (xyz) in the 
shaded stripe S{ti,t2) carry the topological term iOqi^^y^y 



various deficiencies, as was shown by very accurate Monte 
Carlo simulations of the step scaling function introduced 
in [2^, the 2-d 0(3) model with a topological lattice ac- 
tion still has the correct quantum continuum limit ^2^. 
Furthermore, although there are even zero-action dislo- 
cations, Xt was found to have only a logarithmic and not 
a power-law divergence. The continuum result for the 
step scaling function is known analytically thanks to an 
ingenious use of the thermodynamic Bethe ansatz (23| . 
Remarkably, the step scaling function is known analyti- 
cally also at 6 = TT [2J|. At intermediate values of 9, on 
the other hand, the 2-d 0(3) model is expected not to 
be integrable. Still, by expanding around — n, some 
interesting analytic results have been obtained even in 
that regime [25| . 

As we will see, using very accurate Monte Carlo simu- 
lations, we confirm the analytic results for the step scal- 
ing function at 9 = t: with better than per mille level 
accuracy. This requires good control of lattice artifacts. 
In fact, in the 2-d 0(3) model (at 9 — 0) the behavior 
of lattice artifacts, which is apparently linear in the lat- 
tice spacing a, has been puzzling for many years 26l, |27|. 
Only recently, the puzzle has been resolved by a careful 
analysis in the framework of Symanzik's improvement 
program [ll|, [2^ . It turned out that the apparent linear 
behavior is mimicked by the expected quadratic behavior 
modified by large logarithmic corrections. Interestingly, 
on some lattices the topological action approaches the 
continuum limit of the step scaling function from below, 
while the standard action approaches it from above. Here 
we combine both actions in such a way that cut-off ef- 
fects are extremely reduced to at most a few per mille. 
Using that action as well as the standard action allows 
us to extrapolate reliably to the continuum limit. 

Let us consider the 0(3) model on a triangulated 
square lattice as illustrated in Figure 1, with a 3- 
component unit-vector G S'^ attached to each lattice 
site X. We choose periodic boundary conditions in the 
short spatial direction of extent L and open boundary 
conditions in the long Euclidean time direction. The ac- 
tion is defined on nearest-neighbor bonds (xy) (but not 
on the plaquette diagonals), and is given by 



for Cx ■ Cy > cos (5 and s{ex,ey) = 00 otherwise. This 
action eliminates field configurations for which the angle 
between neighboring spins exceeds 6 [2l|. For 6 = tt the 
constraint becomes irrelevant and the action reduces to 
the standard action, while for 17 = 00 it reduces to the 
topological action of [2^. Besides the standard action, 
we will use an action with an optimized constraint angle, 
cos (5 = —0.345, which has very small cut-off effects 29]. 
Let us also define the geometric topological charge den- 



sity q{xyz) e 



2 ' 2 



associated with a triangle (xyz), 



Rexp{2Triq(^xyz)) = 'i- + Ex ■ ey + Cy ■ + ■ Cx 

+ iex-{eyxez), R>0, (2) 

with 4:TTq(^xyz) being the oriented area of the spherical 
triangle on defined by the three unit- vectors et^, Cy, 
and Cz- On a periodic lattice, the topological charge 
Q = '^{xyz) 1{xyz) would be an integer in the second ho- 
motopy group 112 [S*^] — ^- Here we work with open 
boundary conditions in Euclidean time and we sum the 
topological charge density only over the stripe S{ti,t2) 
of shaded plaquettes between ti and ^2 illustrated in 
Figure 1. This yields the non- integer valued quantity 
Q{ti,t2) = J2{xyz)es{tut2)'ii^y^)- I" *° determine 

the massgap, we introduce the operator E{t) — J2x 
where the sum extends over all points x = ixi,t) in a 
time-slice. We now define the 2-point function 
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'S'[e| = ^ s(e*c,ej,), s{ex,ey) = ^{1 



X (1) 
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X exp(-S'[^ + i9Q{ti,t2)) - cxp(-m(6i, L){t2 - h)), 
Ziti,t2;9)=ll [dex eM~S[e\+i9Q{ti,t2)), (3) 

which decays exponentially with the 9- and L-dependent 
massgap m{9, L) at large Euclidean time separations. 
The massgap has been determined with high accuracy 
from numerical simulations using the meron-cluster al- 
gorithm as well as a variant of a method developed 
by Hasenbusch [lo| , which was further improved in j28} . 
It is straightforward to include the topological term in 
this method. Numerically, we obtain C{ti,t2;9) as the 
ratio ofC{ti,t2;9)Z{ti,t2;9)/Z{0) and Z{ti,t2;9)/Z{0), 
where Z{0) — Z(ti, t2] 0), which is independent of <i and 
^2. It turns out that the resulting complex action prob- 
lem for ^ is mild for the moderate spatial volumes 
that are relevant in this study. 

Following 22|, we set — Lm{9,L) and define the 
step 2 scaling function T,{9,2,uo,a/ L) = 2Lm{9, 2L). 
Remarkably, both for e = fill and for 6* = tt [2^, 
the continuum limit a{9,2,uo) — S(0, 2, itg, a/L — > 0) 
has been determined analytically. One particular value 
is a{9 = 0,2,uo = 1.0595) = 1.2612103 [2^. Figure 2 
shows the cut-off effects of I](0, 2, uq, a/L) for the stan- 
dard action, the topological action, and the optimized 
constraint action. The constraint has been optimized 
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FIG. 2. Cut-off dependence of the step scaling function 
E(S, 2,uo = 1.0595, a/L) for three different lattice actions: 
the standard action, the topological lattice action of J2dl l. and 
the optimized constraint action with cos 5 = —0.345, as well 
as for three different values of 6 = ( top ) and 6 = (bot- 
tom). The lines are fits based on eq.^. The horizontal lines 
represent the analytic continuum results for 9 = J2^] and 
9 = n fKil . and the fitted continuum value for 9 = ^. 



to COS 5 = —0.345 by demanding that E(0, 2, uo, a/i = 
= a(0,2,uo) for uq = 1.0595. For L/a > 6, the 
remaining cut-off effects are then less than a per mille. 
The data are fitted to Synianzik's effective theory, which 
predicts ^ 



S(0,2,wo,a/i) ^ 
L2 



(t(0,2,wo) 



[Blog\L/a) + C\og'{L/a) 



■(4) 



They are in excellent agreement with the analytic pre- 
diction to four significant digits accuracy. 

Figure 2 also shows corresponding results a.t 9 = it 
and 7r/2, both for the standard action and for the op- 
timized constraint action with cos (5 — —0.345. At 
9 = TT the analytic result for the step scaling function 
is cr(7r,2,uo = 1.0595) = 1.231064 f22|. Remarkably, al- 
though the constraint was optimized for 9 — 0, the cut-off 



effects are still at most a few per mille at 6* = tt. Since 
exp(i7r(5) = (—1)*^ = exp(— jttQ), ior 9 = tt the topologi- 
cal term does not explicitly break any symmetries of the 
9 = theory. Consequently, the corresponding analysis 
of the cut-off effects |28j] , which underlies the fits in Fig- 
ure 2, still applies. Again, we find excellent agreement 
with the analytic result to four digits accuracy, which in- 
directly confirms the corresponding exact S-matrix 12 1 
that underlies the analytic calculation. Agreement of the 
same quality is obtained for other values of uq . 

Finally, let us consider the case ^ = § , which is again il- 
lustrated in figure 2. For 9 ^ or n there are no analytic 
results. The topological term then breaks both parity 
and charge conjugation and hence the detailed analysis of 
the cut-off effects in no longer applies. Still, remark- 
ably the optimized constraint action with cos 5 = —0.345 
has no observable cut-off effects, and extrapolates to the 
same continuum value ct(^,2, uq = 1.0595) = 1.24733(4) 
as the standard action. This value of the step scaling 
function differs from the values at 6* = and 9 = tt hy 
three significant digits. This shows that the continuum 
limit a.t 9 = ^ indeed represents a different theory. The 
same is true for other intermediate values < 9 < n. 
Hence, 9 is indeed a relevant physical parameter that 
does not renormalize to or tt non-perturbatively, as one 
might have expected due to the presence of dislocations. 

If xt is logarithmically ultra-violet divergent in the con- 
tinuum limit, which may not be the case for other defini- 
tions of the topological charge [s^ , the difference between 
the energy densities of different 0- vacua diverges as well. 
This is a peculiarity of the 2-d 0(3) model, which should 
not extend to CP{N—1) models with N > 2. Despite the 
divergence of the vacuum energy density, the massgap 
m{9, L) is completely well-behaved, and a proper non- 
trivial continuum field theory can be defined for any value 
oi9. It will be interesting to further investigate these the- 
ories. At large volume L and low energies, for = tt one 
expects the 0(3) symmetry to dynamically enhance to 
the 0(4) = SU{2)l x SU{2)r symmetry of the WZNW 
model. In the infinite volume limit, the triplet state 
should then become massless, i.e. m{9 = tt,L oo) 0, 
and degenerate with an 0(3) singlet state. For small val- 
ues of 9, the singlet state is a scattering state with an 
energy 2171(9, L — >■ oo). For example, it would be inter- 
esting to figure out whether there is a critical value 9c 
at which the singlet becomes a bound state 25|. Since 
9 does not get renormalized, despite the fact that the 
energy density difference between different 0-vacua may 
be logarithmically ultra-violet divergent, the massgap is 
a finite physical quantity, whose value at = tt is even 
known analytically; e.g. at Lm{0, L) = 1.0595 one ob- 
tains Lni{-K,L) = 1.048175 As shown in Figure 3, 
the numerical data for L = 24a, which extend to ar- 
bitrary values of 9, agree with this analytic prediction 
below the permille level. There is a remaining tiny cut- 
off effect, which diminishes only for volumes as large as 
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FIG. 3. The d-dependent massgap Lm{6,L) at Lm{Q,L) — 
1.0595 using the optimized constraint action for L = 24a, 
compared to the analytic result at = n jfej / (cross). 



L = 128a. The cut-ofF effects of constraint actions will 
be discussed in detail in forthcoming publications [29]. 

Since dislocations do not prevent a non-trivial contin- 
uum limit in the 2-d 0(3) model with 6* 7^ 0, the same is 
expected for 2-d CP{N — 1) models. Using an unconven- 
tional D-theory regularization by SU{N) quantum spin 
ladders, <CP(N — 1) models have been simulated success- 
fully at 9 — TT [3l|. More recently, a worm algorithm 
has been constructed for CP{N — 1) models [33], which, 
however, suffers from a sign problem at 9 ^ 0. Still, this 
algorithm or a variant of Hasenbusch's method may be 
sufficiently powerful to simulate CP{N — 1) models for 
arbitrary values of 6, and thus generalize the results we 
have obtained here. Furthermore, it is obvious to ask 
whether 6*- vacua effects in 4-d Yang-Mills theories can be 
addressed with similar methods, using a geometric def- 
inition of the lattice topological charge 33-i3(| . In that 
case, a potential dislocation problem arises for SU{2) and 
SU{3) i37|, but can be cured by using an improved lattice 
action |38|. Based on the results obtained here, we ex- 
pect that dislocations are harmless also in 4-d Yang-Mills 
theories, and that 0-vacuum effects can be simulated at 
least in moderate spatial volumes. 

We dedicate this paper to P. Hasenfratz on the occa- 
sion of his 65th birthday. Over many years, we have ben- 
efitted tremendously from his deep insights into strongly 
interacting field theories, which he generously shared 
with us. We are indebted to J. Balog, M. Liischer, and 
P. Weisz for illuminating discussions, and to J. Balog 
for providing exact values of the step scaling function 
at = TT prior to publication. This work has been 
supported by the Regione Lombardia and CILEA Con- 
sortium through a LISA 2011 grant, as well as by the 
Schweizerischer Nationalfonds (SNF). 
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